Mon.  Not.  R.  Astron.  Soc.  345,  1255-1265  (2003) 


Thermal  and  non-thermal  plasmas  in  the  galaxy  cluster  3C  129 

H.  Krawczynski,1*  D.  E.  Harris,2  R.  Grossman,3  W.  Lane,4  N.  Kassim4 
and  A.  G.  Willis5 

1  Washington  University  in  St.  Louis.  Phys.  Department,  I  Brookings  Drive,  Campus  Box  1105,  St.  Louis,  MO  63130,  USA 
2 Smithsonian  Astrophysical  Observatory,  60  Garden  Str.,  Cambridge,  MA  02138,  USA 
3  Department  of  Physics  and  Astronomy.  Barnard  College,  Columbia  University,  New  York,  NY  10027,  USA 
4 Code  7213,  Remote  Sensing  Division,  Naval  Research  Laboratory,  Washington,  DC  20375-5351,  USA 

5  National  Research  Council  of  Canada,  Herzberg  Institute  of  Astrophysics,  Dominion  Radio  Astrophysical  Obserx’atory,  PO  Box  248, 
Penticton,  BC,  Canada  V2A  6K3 


Accepted  2003  July  23.  Received  2003  July  22;  in  original  form  2003  January  28 


ABSTRACT 

We  describe  new  Chandra  spectroscopy  data  of  the  cluster  which  harbours  the  prototypical 
‘head-tail’  radio  galaxy  3C  129  and  the  weaker  radio  galaxy  3C  129.1.  We  combined  the 
Chandra  data  with  Very  Large  Array  (VLA)  radio  data  taken  at  0.33,  5  and  8  GHz  (archival 
data)  and  1.4  GHz  (new  data).  We  also  obtained  new  H  I  observations  at  the  Dominion  Ra¬ 
dio  Astrophysical  Observatory  (DRAO)  to  measure  the  neutral  hydrogen  column  density  in 
the  direction  of  the  cluster  with  arcmin  angular  resolution.  The  Chandra  observation  reveals 
extended  X-ray  emission  from  the  radio  galaxy  3C  129.1  with  a  total  luminosity  of  1.5  x 
1041  erg  s'1.  The  X-ray  excess  is  resolved  into  an  extended  central  source  of  ~2  arcsec 
(1  kpc)  diameter  and  several  point  sources  with  individual  luminosities  up  to  2.1  x  1040  erg 
s' 1 .  There  is  no  evidence  for  a  correlation  between  the  3C  1 29. 1  X-ray  and  radio  morphology. 
In  the  case  of  the  radio  galaxy  3C  129,  the  Chandra  observation  shows,  in  addition  to  core 
and  jet  X-ray  emission  reported  in  an  earlier  paper,  some  evidence  for  extended,  diffuse  X-ray 
emission  from  a  region  east  of  the  radio  core.  The  12  x  36  arcsec2  (6  kpc  x  17  kpc)  region  lies 
‘in  front’  of  the  radio  core,  in  the  same  direction  into  which  the  radio  galaxy  is  moving.  We  use 
the  radio  and  X-ray  data  to  study  in  detail  the  pressure  balance  between  the  non-thermal  radio 
plasma  and  the  thermal  intracluster  medium  (ICM)  along  the  tail  of  3C  129  which  extends  over 
15  arcmin  (427  kpc).  Depending  on  the  assumed  lower  energy  cut-off  of  the  electron  energy 
spectrum,  the  minimum  pressure  of  the  radio  plasma  lies  a  factor  of  between  10  and  40  below 
the  ICM  pressure  for  a  large  part  of  the  tail.  We  discuss  several  possibilities  to  explain  the 
apparent  pressure  mismatch. 

Key  words:  galaxies:  active  -  galaxies:  clusters:  individual:  3C  129  -  galaxies:  jets  -  radio 
continuum:  galaxies  -  X-rays:  galaxies  -  X-rays:  galaxies  -  clusters. 


1  INTRODUCTION 

Observations  of  radio  galaxies  within  galaxy  clusters 

Tailed  radio  galaxies  within  galaxy  clusters  have  drawn  intensive 
study  for  several  decades  on  both  theoretical  and  observational 
grounds.  Radio  and  X-ray  data  complement  each  other  in  that  the  ra¬ 
dio  data  show  high  resolution  images  of  the  jets  and  more  extended 
synchrotron  features  yielding  information  on  strength  and  orienta¬ 
tion  of  the  magnetic  field  within  the  radio  plasma  and  between  the 
radio  galaxy  and  the  observer,  while  the  X-ray  data  gives  information 
about  the  temperature,  chemical  composition,  density,  and  pressure 
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of  the  intracluster  medium  (ICM).  Thus,  the  jets  propagate  into  a 
well-defined  medium,  which  greatly  facilitates  the  interpretation  of 
the  radio  images.  It  becomes  possible  to  study,  for  example,  the  pres¬ 
sure  balance  between  the  radio  plasma  and  ICM.  and  thus  to  search 
for  hidden  pressure  components  in  the  radio  plasma.  Various  authors 
have  studied  the  radio  plasma/ICM  pressure  balance  for  FR I  and  FR 
II  radio  galaxies,  in  most  cases  finding  the  minimum  synchrotron 
pressure  lower  than  the  ICM  pressure  by  approximately  one  order  of 
magnitude  (e.g.  Feretti  et  al.  1990;  Hardcastle  &  Worral  2000,  refer¬ 
ences  therein).  Furthermore,  the  long  tails  of  radio  galaxies  carry  the 
imprint  of  relative  motion  between  the  non-thermal  plasma  and  the 
ambient  hot  gas.  In  the  parlance  of  the  field,  they  reflect  the  weather 
conditions  in  the  ICM.  The  potential  of  such  observations  to  reveal 
details  of  cluster  mergers  such  as  subsonic  and  transonic  bulk  flows, 
shocks,  and  turbulence  has  been  investigated  on  theoretical  grounds 
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Table  1.  3C  129  Chandra  data  sets. 


Description 

Energy 

Integration  time 

Ang.  res.  (FWHM) 

Reference 

Chandra  ACIS  I 

0.5-8  keV 

9.3  ksec 

0.84  arcsec  (aimpoint) 

Krawczynski  (2002) 

Chandra  ACIS  S 

0.5-8  keV 

28.3  ksec 

0.84  arcsec  (aimpoint) 

Krawczynski  (2002) 

by  a  number  of  authors  (Loken  et  al.  1995;  Burns  et  al.  2002).  On 
the  observational  side,  Gomez  et  al.  (1997)  studied  a  large  number 
of  wide-angle  tail  (WAT)  radio  galaxies  inside  Abell  clusters.  They 
found  a  strong  correlation  between  the  orientation  of  the  radio  tails 
and  the  elongations  of  the  ICM  surface  brightness  profiles,  and  in¬ 
terpreted  it  as  evidence  for  bulk  plasma  motion  that  follows  cluster 
mergers.  Bending  of  radio  tails  has  often  been  explained  by  invok¬ 
ing  motion  of  the  radio  galaxy  in  combination  with  ‘stormy  ICM 
weather’  (e.g.  Feretti  et  al.  1985;  Sakelliou,  Merrifield  &  McHardy 
1996)  but  this  interpretation  does  not  work  for  all  bent  radio  tails 
(e.g.  Eilek  et  al.  1984). 

In  addition  to  facilitating  studies  of  radio  galaxies  and  the  ICM 
themselves,  observations  of  radio  galaxies  in  galaxy  clusters  make  it 
possible  to  study  the  interaction  of  the  radio  galaxies  with  the  ICM 
and  vice  versa.  The  radio  plasma  can  interact  with  the  ICM  and  form 
so-called  cavities  devoid  of  ICM  that  show  up  as  depressions  in  the 
X-ray  surface  brightness  maps.  Cavities  have  been  observed  in  the 
Perseus  cluster  Bohringer  et  al.  (1993),  in  the  Cygnus  A  cluster 
by  Carilli,  Perley  &  Harris  (1994);  more  recently,  Chandra  obser¬ 
vations  revealed  cavities  in  the  clusters  Hydra  A,  A2597,  A2052 
and  A4059  (McNamara  et  al.  2000;  2001;  Blanton,  Sarazin  &  Irwin 
2001b;  Heinz  et  al.  2002).  It  is  not  yet  clear  whether  these  cavities 
are  merely  curiosities  or,  whether  they  influence  the  ICM  in  an  im¬ 
portant  way  by  supplying  a  dynamically  important  amount  of  heat 
and/or  magnetic  field  to  the  ICM  (Briiggen  &  Kaiser  2002;  Quilis, 
Bower  &  Balogh  2001;  Reynolds,  Heinz  &  Begelman  2002;  Fabian 
et  al.  2002).  The  relative  motion  of  the  ICM  and  the  Interstellar 
Medium  (ISM)  of  the  radio  galaxy’s  host  galaxy  can  result  in  a 
modification  of  the  ISM  and  in  a  bow  shock  forming  ‘in  front"  of 
the  radio  galaxy  (e.g.  Stevens,  Acreman  &  Ponman  1999;  Toniazzo 
&  Schindler  2001,  and  references  therein).  Although  these  interac¬ 
tions  can  also  be  observed  in  normal  galaxies,  the  radio  maps  of 
radio  galaxies  can  be  used  to  interpret  the  data  better  by  constrain¬ 
ing  the  motion  of  the  galaxies  in  the  plane  of  the  sky  as  well  as  the 
ICM  magnetic  field. 

The  cluster  3C  129  and  outline  of  the  paper 

In  this  paper  we  present  Chandra  data  on  the  nearby  galaxy  clus¬ 
ter  3C  129  (z  —  0.0223).  Owing  to  its  low  Galactic  latitude 
(/  =  16075,  b  —  073)  it  has  not  been  studied  extensively  at  optical 
wavelengths.  In  addition  to  velocity  measurements  of  3C  129  and 
3C  129.1  (Spinrad  1975),  we  are  aware  of  only  one  other  velocity- 
verified  cluster  galaxy:  WEIN  048  (Nilsson  et  al.  2000).  Optically, 
the  galaxy  3C  129  has  been  described  as  a  weak  elliptical  galaxy 
without  peculiar  features  (Colina  &  Perez-Fournon  1990).  X-ray  ob¬ 
servations  with  the  Einstein,  EXOSAT,  and  ROSAT  satellites  (Edge 
&  Stewart  1991;  Leahy  &  Yin  2000)  yielded  a  cluster  luminosity  of 
2.7  x  1044  erg  s-1  in  the  0.2-10  keV  energy  band,  a  central  gas  den¬ 
sity  of  6  x  10-3cm-3  and  an  ICM  mass  of  3.6  x  1013  Mq.  Based  on 
Einstein  IPC  data  Leahy  &  Yin  (2000)  found  evidence  for  a  cooling 
flow  with  M  ^  85  Mq  yr^1.  However,  reanalyzing  the  same  data, 
Taylor  et  al.  (2001)  derived  an  upper  limit  of  M  <  1.2  Mq  yr_1. 
Based  on  the  Chandra  imaging  spectroscopy  data  we  performed  a 


spectroscopic  deprojection  analysis,  taking  into  account  the  asym¬ 
metry  of  3C  129's  X-ray  surface  brightness  by  using  a  set  of  elliptical 
ICM  shells  (Krawczynski  2002,  called  Paper  I  in  the  following).  This 
analysis  did  not  show  a  decrease  of  the  ICM  temperature  toward  the 
cluster  centre,  as  expected  for  a  cooling  flow. 

The  cluster  contains  two  radio  galaxies.  The  radio  source  3C  129 
is  a  well  studied  galaxy  in  the  radio  wavelength  range  as  a  prototype 
head-tailed  (HT)  galaxy  with  a  strongly  curved,  15  arcmin  long 
(426  kpc),  two-armed  tail  (Miley  1973;  van  Breugel  &  lagers  1982; 
Kassim  et  al.  1993;  Feretti  et  al.  1998).  In  an  earlier  paper  (Harris, 
Krawczynski  &  Taylor  2002,  hereafter  Paper  II)  we  presented  ev¬ 
idence  for  X-ray  emission  from  the  inner  northern  jet  of  3C  129. 
Near  the  projected  centre  of  the  X-ray  galaxy  cluster  is  the  weaker 
and  smaller  radio  source  3C  129.1 .  This  galaxy  also  has  dual  radio 
jets  which  extend  over  —2  arcmin  (57  kpc)  (Downes  1984;  Jagers 
1987;  Kassim  et  al.  1993). 

In  this  paper  we  describe  X-ray  and  radio  observations  of  this 
cluster.  We  complement  the  Chandra  X-ray  data  with  archival  0.33, 
5  and  8  GHz  and  new  1.4-GHz  Very  Large  Array  (VLA)1  radio 
data.  We  use  the  radio  and  X-ray  data  to  search  for  evidence  of 
interactions  between  the  radio  plasma  and  the  ICM  and  to  study 
the  pressure  balance  between  the  two  plasmas.  Since  the  cluster 
lies  at  low  Galactic  latitude  we  expect  substantial  variations  of  the 
neutral  hydrogen  column  density  Ah  over  the  field  of  view  of  the 
Chandra  observations.  We  therefore  obtained  H I  observations  with 
the  DRAO  telescope.  Tables  1  and  2  summarize  the  data  sets  used 
in  our  analysis. 

The  remainder  of  the  paper  is  structured  as  follows.  We  introduce 
observations  and  data  sets  in  Section  2.  Results  of  radio  observations 
are  presented  in  Section  3,  while  in  Section  4  we  discuss  radial 
profiles  of  the  ICM  surface  brightness  distribution  together  with  the 
geometry  of  the  ICM.  We  describe  the  search  of  the  X-ray  data 
for  signatures  of  radio  galaxy /ICM  interactions  in  Section  5,  and 
examine  the  pressure  balance  between  the  plasma  of  3C  129  and  the 
ICM  in  Section  6.  Finally,  we  summarize  our  results  in  Section  7. 

We  use  Ho  =  ho  100  km  s_1  Mpc-1  with  ho  =  0.65  throughout, 
setting  the  cluster  at  a  distance  of  100  Mpc;  1  arcsec  corresponds 
to  474  pc.  We  quote  all  errors  on  90  per  cent  confidence  level.  All 
coordinates  are  given  for  equinox  2000.  Distances  from  the  radio 
core  of  3C  129  are  given  as  point-to-point  distances,  and  not  as 
distances  along  the  radio  tail. 

2  DATA  SETS  AND  DATA  ANALYSIS 

2.1  H  i  observations  with  the  Dominion  Radio 
Astrophysical  Observatory 

Since  the  3C  129  area  lies  close  to  the  Galactic  plane,  it  is  included 
in  the  21 -cm  line  and  continuum  observations  planned  as  part  of  the 

1  The  National  Radio  Astronomy  Observatory  is  a  facility  of  the  National 
Science  Foundation  operated  under  cooperative  agreement  by  Associated 
Universities,  Inc. 
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Table  2.  Radio  data  sets  used  in  this  paper. 


Frequency 

Instrument 

New  data 

Ang.  res.  (FWHM) 

Noise  RMS  [mJy  beam  1 

1  DR 

Reference 

Hi 

DRAO 

Y 

1  arcmin 

- 

- 

this  work 

330  MHz 

VLA 

N 

8  arcsec 

0.75 

400 

Lane  et  al.  (2002) 

1.4  GHz 

VLA 

Y 

16  arcsec 

0.1 

3000 

this  work 

4.7  GHz  (3C  129) 

VLA 

N 

1.8  arcsec 

0.072 

545 

Taylor  et  al.  (2001) 

8  GHz  (3C  129) 

VLA 

N 

1.8  arcsec 

0.04 

925 

Taylor  et  al.  (2001) 

8  GHz  (3C  129) 

VLA 

N 

0.83  arcsec 

0.03 

1287 

Harris  et  al.  (2002) 

8  GHz  (3C  129.  1) 

VLA 

N 

0.8  arcsec 

0.012 

225 

Taylor,  private  communication 

second  phase  of  the  Canadian  Galactic  Plane  Survey  (CGPS;  see 
Taylor  1999,  for  an  overview  of  its  goals).  The  DRAO  Synthesis 
Telescope  H  I  observations  of  this  part  of  the  sky  were  completed 
early  in  the  second  phase  of  the  survey.  They  were  reduced  follow¬ 
ing  the  standard  H  I  reduction  procedure  for  the  CGPS,  described 
by  Higgs  (1999).  This  work  includes  the  insertion  of  ‘short  spac¬ 
ing  data’  collected  with  the  DRAO  26-m  telescope  and  corrected 
for  stray  radiation  effects.  So  the  H  I  data  are  fully  sampled  out 
to  an  interferometer  spacing  of  approximately  2900  wavelengths. 
The  resulting  images  have  a  resolution  of  1  arcmin  full  width  at 
half-maximum  (FWHM)  in  right  ascension  and  1  arcmin  cosec(S) 
in  declination  5.  We  used  a  method  similar  to  that  described  by 
Normandeau  (1999)  to  determine  the  optical  depth  in  the  direction 
toward  the  3C  129  cluster.  Briefly,  H  I  spectra  in  the  direction  of 
the  heads  of  the  radio  galaxies  3C  129  and  3C  129.1,  which  show 
absorption,  were  compared  with  spectra  toward  nearby  ‘off-source’ 
positions.  The  comparisons  indicate  that  the  optical  depth  along  the 
line  of  sight  to  these  two  strong  sources  averages  to  less  than  1. 
We  assume  that  this  optically  thin  condition  holds  anywhere  in  the 
direction  toward  the  3C  129  cluster,  so  we  can  derive  the  H I  column 
density  from  the  expression: 

Ah  =  1.823  x  1018  x  Y.iTb(vi)Avi.  (1) 

Here,  (VH  is  given  in  atom  cm-2,  the  sum  runs  over  all  velocities 
Vj  which  contain  Galactic  hydrogen,  Tb  is  the  radio  brightness  tem¬ 
perature  in  degrees  Kelvin,  and  A  u,-  is  the  width  of  the  ith  velocity 
interval  in  km  s-1. 

2.2  Radio  observations  with  the  Very  Large  Array 

The  330-MHz  VLA  data  were  acquired  during  175  snapshot  obser¬ 
vations  made  over  several  years  in  the  late  1990s.  The  data  were 
taken  in  A,  B  and  C  configurations.  The  final  image  has  ^8  arcsec 
resolution  (FWHM).  The  root  mean  square  (RMS)  noise  in  the  map 
is  ~0.75  mJy  beam-1,  and  the  dynamic  range  (DR,  the  signal-to- 
noise  ratio  of  the  peak  brightness)  is  ~400.  A  detailed  description 
of  these  observations  has  been  given  by  Lane  et  al.  (2002). 

Data  for  the  1400-MHz  image  were  obtained  at  the  VLA  on  2001 
August  16  (4  h,  C  configuration),  and  2001  December  6  (1.5  h,  DnC 
Configuration).  Short  observations  of  3C  48  were  used  for  amplitude 
calibration,  and  the  data  were  combined  and  mapped  with  the  pack¬ 
age  Alps  (Astronomical  Imaging  Processing  System)  using  standard 
imaging  techniques.  The  image  has  a  FWHM  resolution  of  «U8  x 
14  arcsec2  at  a  position  angle  of  —15°,  an  RMS  noise  of  100  pJy 
beam-1,  and  a  DR  ~3000.  The  primary  aim  of  the  VLA  observa¬ 
tions  in  the  L  band  was  to  derive  accurate  0.33-1.4  GHz  spectral 
indices.  The  330  MHz  data  have  a  minimum  UV-spacing  of  40A, 
and  a  maximum  of  40.65  kA  (k  stands  for  kilo);  for  the  1400-MHz 
data,  the  minimum  spacing  is  180A  and  the  maximum  is  16.5  kA.  To 
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make  the  spectral  index  maps,  both  data  sets  were  limited  to  a  maxi¬ 
mum  UV  of  15  kA  and  a  Gaussian  taper  was  applied  during  imaging 
to  match  the  measured  beam  sizes  as  closely  as  possible.  A  circular 
restoring  beam  of  18  x  18  arcmin2  was  used  in  the  final  maps.  We 
estimated  the  errors  on  the  spectral  indices  by  adding/subtracting 
5ff  to  fluxes  at  one  frequency  and  subtracting/adding  5<r  from  the 
other.  This  method  gives  conservative  5 a  errors  of  0.01  near  the 
radio  core  and  0.25  at  the  end  of  the  radio  tails. 

Finally  we  used  several  VLA  maps  of  3C  129  and  3C  129.1  with 
various  effective  angular  resolutions  that  G.  Taylor  kindly  provided. 
The  pressure  determination  of  the  3C  129  radio  plasma  uses  a  5— 
8  GHz  spectral  index  map  extracted  from  two  data  sets  with  identical 
angular  resolutions  of  1.8  arcsec  (Taylor  et  al.  2001).  The  width  of 
the  3C  129  tail  was  determined  at  distances  from  0  to  20  arcsec, 
0.3  to  2.5  arcmin,  and  2.5  to  15  arcmin  from  the  radio  core  from 
the  8-GHz,  0.83-arcsec  angular  resolution  map  from  Harris  et  al. 
(2002),  the  8-GHz,  1.8-arcsec  angular  resolution  map  from  Taylor 
et  al.  (2001),  and  the  0.33-GHz,  8-arcsec  angular  resolution  map 
from  Lane  et  al.  (2002),  respectively. 

A  high  angular  resolution  map  of  3C  129.1  (0.81"  x  0.76" 
FWHM)  has  been  used  for  the  radio  contours  shown  later  (in  Fig.  5). 


2.3  Chandra  observations 

Our  data  consist  of  two  pointings:  a  3 1 .4-ks  exposure  taken  on  2000 
December  9  with  the  ACIS  S,  6-chip  array  (aimpoint  close  to  3C 
129),  and  a  9.6-ks  exposure  with  the  standard  ACIS  I  CCD  on  2001 
January  9  (aimpoint  close  to  3C  129.1).  The  analysis  uses  the  ciao 
2.2  software.  Filtering  for  major  background  flares  reduces  the  ACIS 
S  data  set  to  28.3  ks  and  the  ACIS  I  data  set  to  9.3  ks.  We  produce 
flat-fielded  X-ray  brightness  maps  using  a  soft  (0.5-1  keV),  an  in¬ 
termediate  (1-2  keV)  and  a  hard  (2-5  keV)  energy  band.  The  bright¬ 
ness  maps  were  corrected  for  instrumental  parameters  such  as  the 
effective  area  etc,  by  dividing  by  the  appropriate  exposure  maps. 
The  flat-fielded  maps  give  the  flux  of  0.5-5  keV  photons  cm-2  s 
arcsec-2.  We  combine  the  flat-fielded  maps  of  both  observations  by 
weighting  each  one  with  the  corresponding  integration  time.  Further 
details  concerning  these  data  can  be  found  in  Paper  I. 

Some  X-ray  features  have  only  modest  statistical  significance. 
Using  the  number  of  counts  in  the  source  ‘ON  region’  A0n.  and 
the  number  of  counts  in  the  comparison  ‘OFF  region"  N off.  we 
compute  the  statistical  significances  based  on  the  method  of  Li  & 
Ma  (1983).  The  method  is  based  on  Poisson  probability  distributions 
and  correctly  accounts  for  the  statistical  uncertainty  of  the  counts  in 
both  the  source  and  comparison  regions.  The  statistical  significance 
of  an  excess  is  given  in  standard  deviations  by 

S(Aon>  A0fG  ?)  =  \/— 21n(A),  (2) 
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and  f  given  by 
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where  A£20n  and  A£20ff  are  the  solid  angle  areas  of  the  source 
and  background  regions,  respectively,  and  A(£2)  is  the  the  effective 
area  of  the  detector.  We  use  the  exposure  map  produced  for  1.7- 
keV  photons  as  A(£2)  and  integrate  over  this  map  with  the  ciao 
tool  dmstat.  Although  the  method  is  based  on  Poisson  probability 
distributions,  the  expression  for  the  statistical  significance  is  only 
valid  for  A0 N  >  10  and  A  off  >  10  (Alexandreas  et  al.  1992).  If 
this  condition  is  not  met,  we  use  the  equation  of  Alexandreas  et  al. 
(1992), 
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to  compute  the  chance  probability  Pc  that  A0 N  counts  are  observed 
by  chance,  taking  into  account  the  statistical  uncertainty  of  the  ‘OFF’ 
measurement.  We  state  the  significances  without  taking  into  account 
‘trial  factors’  owing  to  the  search  for  features  over  an  extended  solid 
angle  region. 


3  RADIO  MAPS 

The  Ah  map  derived  from  the  H  I  observations,  smoothed  with  a 
Gaussian  of  2.45  arcmin  FWHM  is  shown  in  Fig.  1.  Although  the 
AH  values  change  by  25  per  cent  over  the  1725  x  2? 5  field  of  view. 
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Figure  1.  Neutral  hydrogen  column  density  in  the  direction  of  the  cluster 
3C  129  derived  from  the  DRAO  H I  data,  smoothed  with  a  Gaussian  of  2.45 
arcmin  FWHM.  Darker  regions  show  higher  Ah  column  densities.  Rather 
atypically,  the  column  densities  are  a  bit  lower  along  the  Galactic  plane  than 
at  slightly  lower  and  higher  latitudes.  The  ten  contour  levels  are  linearly 
spaced  from  8  x  1021  to  1022  cm-2.  Compared  with  nearby  fields  in  the 
Galactic  plane  (see  e.g.  Leahy  &  Roger  1991,  figs  7  and  8),  only  modest 
gradients  of  the  Ah  column  density  are  observed. 
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the  Chandra  observations  lie  in  a  region  of  rather  constant  values 
of  between  8.4  x  1021  cm'2  and  9.2  x  1021  cm-2. 

Strasser  &  Taylor  (2002)  studied  H  I  absorption  spectra  of  273 
background  sources  from  the  Canadian  Galactic  Plane  Survey.  In 
their  fig.  2  they  show  that  in  optically  thin  regions  the  apparent  off- 
source  temperatures  cover  a  range  of  about  100  K  (upper  bound) 
to  50  K  (lower  bound).  Assuming  a  mean  value  of  75  K,  these  off- 
source  data  sets  imply  a  maximum  error  in  the  H  I  determination  of 
about  30  per  cent.  Our  results  are  marginally  consistent  with  those 
of  Leahy  &  Yin  (2000)  and  Taylor  et  al.  (2001),  derived  from  fitting 
the  cluster’s  X-ray  energy  spectra:  the  first  authors  report  Ah  = 
(5.7  ±  0.3)  x  1021  cm'2  based  on  ROSAT  and  EXOSAT  data,  and 
the  latter  authors  get  Ah  =  (6.4  ±  0.4)  x  1021  cm'2  from  ASCA 
data. 

The  1 .4-GHz  radio  surface  brightness  distribution  of  the  3C  1 29 . 1 
and  3C  129  region  is  shown  in  Fig.  2,  upper  panel.  The  map  shows 
substantially  more  detail  than  the  previously  published  map  at  this 
frequency  (Jagers  1987).  One  clearly  recognizes  the  peculiar  radio 
feature  (called  ‘radio  crosspiece’  in  the  following)  to  the  northeast  of 
the  radio  head  seen  in  earlier  low-frequency  0.33  GHz  data  (Jagers 
&  de  Grijp  1983;  Lane  et  al.  2002).  Fig.  2,  lower  panel,  shows  the 
0.33-1 .4  GHz  spectral  indices  (a  defined  by  Sv  <xv~a)  derived  from 
the  0.33-GHz  data  of  Lane  et  al.  (2002)  and  our  1.4-GHz  data.  The 
spectrum  is  flattest  at  the  3C  129  radio  core  (acore  =  0.5).  Between 
0  and  6  arcmin  from  the  radio  core  the  spectrum  of  both  radio  tails 
stays  rather  constant  with  spectral  indices  from  0.5  to  0.8;  between 
6  and  1 1  arcmin  from  the  radio  core  it  steepens  to  a  sa  1 ;  and  from 
11  to  15  arcmin  it  varies  from  1  to  1.5.  The  region  of  the  radio 
crosspiece  has  spectral  indices  of  between  1.1  and  1.5.  A  natural 
interpretation  of  the  spectral  index  distribution  in  the  tail  is  that  the 
electron  population  has  an  energy  spectrum  (dA/dL  oc  £'p)  with 
p  as  2.0.  Within  a  distance  of  8  arcmin  from  the  radio  core  the 
radio  spectrum  steepens  from  a  =  0.5  to  1.0,  and  thus  the  electron 
spectrum  steepens  from  p  =  2  to  3,  as  a  consequence  of  synchrotron 
and  inverse  Compton  cooling  of  the  radio  plasma.  Further  cooling 
eventually  results  in  a  >  1,  once  the  electron  spectra  have  cooled 
enough  so  that  the  high-energy  cut-off  of  electrons  is  low  enough 
to  affect  the  frequency  band  being  observed.  Weakening  magnetic 
fields  further  down  the  radio  tail  will  also  contribute  to  this  effect. 

Assuming  no  re-acceleration  of  electrons  in  the  radio  tail  and  us¬ 
ing  an  order  of  magnitude  magnetic  field  strength  of  B  =  5  pG,  we 
estimate  the  age  of  the  electrons  at  8  arcmin  from  the  radio  core 
where  the  spectral  index  of  the  synchrotron  radiation  is  0.5.  For  this 
purpose,  we  assume  an  initial  electron  energy  spectrum  dN /dE 
oc  E~p  with  p  =  2.0,  and  compute  the  0.33-1.4  GHz  spectral  in¬ 
dex  taking  into  account  synchrotron  and  Inverse  Compton  energy 
losses  of  the  electrons  (see  Myers  &  Spangler  1985).  Following 
Jaffe  &  Perola  (1974)  we  assume  that  electrons  are  isotropized  on 
time-scales  short  compared  with  their  radiative  lifetimes.  In  this  ap¬ 
proximation,  the  effect  of  inverse  Compton  cooling  on  the  cosmic 
microwave  background  (CMB)  can  be  taken  into  account  by  us¬ 
ing  an  effective  magnetic  field  B'  =  \J  B2  +  B^mb  with  ^cmb  — 
3.25(1  +  z)2  pG  to  compute  the  energy  losses  of  the  electrons.  The 
model  reproduces  the  observed  spectral  index  a  =  1  for  an  age  of 
95  Myr  of  the  electrons.  This  age  translates  into  a  velocity  of  the 
galaxy  3C  129  of 

ugai  «  2260  sin  (S)-1(/t0/0.65)-1  km  s'1  (6) 

where  6  is  the  angle  between  the  plane  of  motion  of  3C  129  and 
the  line  of  sight  toward  the  observer.  Although  the  equipartition 
magnetic  field  at  8  arcmin  from  the  radio  core  derived  from  the  radio 
data  is  smaller  (between  1 .4  and  3.7  pG  for  lower  frequency  cut-offs 
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Figure  2.  The  upper  panel  shows  the  1.4-GHz  surface  brightness  distribution.  The  peak  surface  brightness  is  0.29  Jy  beam-1  and  the  8  contours  are 
logarithmically  spaced  from  0.001  to  0.256  Jy  beam-1  (beamwidth  16  arcsec  FWHM).  Darker  regions  correspond  to  a  higher  surface  brightness.  The  lower 
panel  shows  the  0.33-1.4  GHz  spectral  indices  (a  from  Sv  oc  v~a)  derived  from  the  0.33-GHz  data  of  Lane  et  al.  (2002)  and  the  1.4-GHz  data  of  the  upper 
panel.  Darker  regions  show  relatively  steeper  energy  spectra  (larger  a)  and  the  7  contour  levels  are  linearly  spaced  from  a  =  0.5  to  a  =  2.  (A a  =  0.25).  The 
small  radio  source  in  the  east  is  3C  129.1  and  the  larger  one  in  the  west  is  3C  129. 


v1  =  330  MHz  and  =  1  MHz,  respectively)  the  radio  plasma  - 
ICM  pressure  balance  argues  for  a  stronger  magnetic  field,  possibly 
up  to  15  pG.  Since  ugai  depends  strongly  on  B ,  a  range  of  allowed 
Ugai  values  constrains  the  magnetic  field  tightly  (see  Section  7). 

Lane  et  al.  (2002)  interpreted  the  radio  crosspiece  as  fossil  non- 
thermal  electrons  re-energized  by  the  bow-shock  of  the  radio  galaxy 
3C  129,  supersonically  traversing  the  ICM.  Assuming  that  the  seg¬ 
ment  of  the  energy  spectrum  of  the  electrons  responsible  for  the 
0.33-  to  1.4-GHz  emission  in  the  radio  crosspiece  has  not  yet  steep¬ 
ened  due  to  radiative  cooling,  the  spectral  index  a  sa  1.3  translates 
into  a  particle  spectral  index  of  p  —  3.6.  Invoking  the  standard  the¬ 
ory  of  diffusive  particle  acceleration  at  shocks  (Drury  1983),  we 
derive  an  estimate  of  the  Mach  number  of  Mgal  =  1.87.  For  an  ICM 
temperature  of  5  keV,  the  ICM  sound  speed  is  1160  km  s-1,  and 
Mgai  =  1-87  corresponds  to  ugai  =  2170  km  s-1  in  fair  agreement 
with  above  estimate.  The  difference  of  the  radial  velocities  of  3C 
129.1  (z  =  0.022265)  and  3C  129  (z  =  0.020814)  is  only  435  km 
s-1,  suggesting  that  the  galaxy’s  velocity  vector  lies  mainly  in  the 
plane  of  the  sky. 

The  spectral  index  values  and  variations  found  across  the  source 
between  data  at  0.33  and  1.4  GHz  are  in  reasonable  agreement  with 


those  found  between  0.61 2  and  1.4GHz  (van  Breugel& lagers  1982; 
Jagers  1987),  and  between  1.4  and  4.1  GHz  (Miley  1973).  Compar¬ 
ing  our  results  with  the  0.074-0.61  GHz  spectral  indices  of  Perley  & 
Erickson  (1979),  we  find  that  the  spectrum  flattens  at  low  frequen¬ 
cies,  while  comparisons  with  the  0.408-2.7  GHz  spectral  indices  of 
Riley  (1973)  and  the  2.7-4.8  GHz  and  4.8-10.6  GHz  spectral  in¬ 
dices  of  Feretti  et  al.  (1998)  show  that  the  energy  spectrum  steepens 
at  higher  frequencies. 


4  THE  GEOMETRY  OF  THE  INTRACLUSTER 
MEDIUM  AND  DISCRETE  SOURCES 

The  upper  panel  of  Fig.  3  shows  the  flat-fielded  data  of  the  two  Chan¬ 
dra  observations,  summed  with  a  weight  proportional  to  the  inte¬ 
gration  time  and  smoothed  with  a  Gaussian  of  18.8  arcsec  FWHM. 
Neither  before  nor  after  smoothing  the  X-ray  surface  brightness 
do  we  detect  obvious  sharp  gradients,  like  those  found  for  a  large 
number  of  more  regular  clusters  (Markevitch,  Vikhlinin  &  Forman 
2002)  which  might  be  produced  by  sloshing  of  the  central  ICM. 
The  cluster  has  a  considerable  ellipticity.  Using  surface  brightness 
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Figure  3.  The  grey-scale  plot  in  the  upper  panel  shows  the  flat-fielded  Chandra  0.5-5  keV  data  (ACTS  I  and  ACIS  S),  smoothed  with  a  Gaussian  of 
18.8  arcsec  FWHM.  Segments  of  ACIS-I  chip  gaps  are  visible  as  white  stripes  outside  the  regions  covered  by  the  ACIS-S  chips.  The  centre  panel  shows  the 
residual  surface  brightness  distribution  after  subtraction  of  the  symmetric  ICM  model;  the  lower  panel  shows  the  residual  surface  brightness  after  subtraction 
of  the  1 8-slice  pie  model.  In  all  three  panels  the  grey-scales  are  linear  and  high  emission  regions  are  darker  than  low  emission  regions;  the  surface  brightness 
in  the  upper  panel  peaks  at  the  location  of  FS  1  with  3.8  x  1 0  6  photon  cm-2  s-1  arcsec-2,  at  the  radio  core  of  3C  129  it  is  1.0  x  10  fl  photon  cm-2  s-1 
arcsec-2;  the  grey  of  the  ‘noise’  at  the  very  end  of  the  3C  129  tail  is  10-7  photon  cm-2  s-1  arcsec-2.  The  centre  and  lower  panels  use  the  same  grey-scales 
with  values  of  2.1  x  10-6  photon  cm-2  s-1  arcsec-2  at  the  location  of  FS  1,  and  10-7  photon  cm-2  s-1  arcsec-2  at  the  noisy  grey  areas  near  the  end  of  the  tail 
of  3C  129,  respectively.  The  radio  contours  have  been  derived  from  the  330-MHz  radio  map  and  indicate  the  location  of  the  radio  galaxy  3C  129. 1  at  the  cluster 
centre  and  3C  129  to  the  east  of  the  cluster  centre.  Contours  are  drawn  at  0.0025,  0.035,  0.08,  0.12,  0.15  and  0.20  Jy  beam-1  (beamwidth  8  arcsec  FWHM). 


isocontours  we  derive  ellipticities  e  =  ( a  —  b)/ a  of  about  0.3,  a  being 
the  semimajor  axis,  and  b  the  semiminor  axis.  The  surface  brightness 
iso-contour  ellipses  are  roughly  aligned  in  the  east-west  direction 
which  depends  slightly  on  the  diameter  of  the  ellipses.  The  surface 


brightness  peaks  at  and  around  the  location  of  two  discrete  sources: 
an  extended  source  near  the  core  of  the  radio  galaxy  3C  129.1  ( a  = 
04h50m06s.6308  and  S  =  45°03'05'.'990),  and  a  point  source,  called 
FS1  in  the  following,  at  2.0  arcmin  toward  the  southwest  of  the 
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Table  3.  X-ray  sources  not  associated  with  known  cluster  galaxies. 


Name 

RA  (Eq.  2000) 
[hh:mm:ss.ss] 

Dec.  (Eq.  2000) 
[dd:mm:ss.ss] 

0.5-5  keV  Fluxa 
[erg  cm-2  s-1] 

Radio  counterpart 

FSl 

FS2 

04:49:58.5 

04:49:54.8 

45:01:44.8 

45:11:31.4 

5.13  x  10-14 
8.85  x  10-14 

0.24  mJy  @4.8  GHz,  0.1 1  mJy  @8  GHz 

no 

“Observed  flux. 


centre.  We  use  the  location  a  =  04h50m02s.9  and  <5  =  45°02'26'.'0 
that  lies  between  these  two  sources  as  the  cluster  centre.  This  lo¬ 
cation  was  determined  by  estimating  centres  of  ellipses  at  various 
surface  brightness  levels. 

Another  strong  point  source,  called  FS2  in  the  following,  is  found 
9.2  arcmin  from  the  cluster  core  toward  its  northwest.  This  second 
source  had  already  been  detected  in  ROSAT  data  by  Leahy  &  Yin 
(2000).  The  location.  X-ray  flux,  and  information  about  radio  coun¬ 
terparts  of  FS1  and  FS2  are  given  in  Table  3.  Lacking  optical  spectra, 
we  can  not  determine  whether  these  two  source  belong  to  the  cluster 
or  not.  The  only  other  known  radio  galaxy  belonging  to  the  cluster, 
WE1N  048  (Nilsson  et  al.  2000),  was  not  in  the  field  of  view  of  the 
two  Chandra  observations. 

We  searched  for  large-scale  irregularities  of  the  ICM  surface 
brightness  by  constructing  models  of  the  undisturbed  ICM  surface 
brightness  and  subtracting  these  model  from  the  surface  brightness 
map.  We  used  two  models  of  the  undisturbed  ICM.  The  first  model  is 
a  ‘standard  King  model'  that  is  symmetric  around  the  cluster  centre: 


E(0  = 


>  +  (r/rcf 


3/3-1 


(7) 


with  the  surface  brightness  E(r),  the  normalization  constant  So, 
the  distance  from  the  cluster  centre  r,  the  cluster  core  radius  rc, 
and  the  exponent  ft.  Using  the  fiducial  value  ft  =  2/3  we  fit  a  core 
radius  of  rc  =  7.1  arcmin  with  a  x2-value  of  317  for  41  degrees 
of  freedom.  Fitting  also  the  /S-parameter  gives  rc  =  7.3  arcmin  and 
ft  =  0.68,  and  does  not  reduce  the  x2  value  substantially:  315  for 
40  degrees  of  freedom.  The  surface  brightness  distributions  after 


subtraction  of  this  King  model  will  show  that  the  ICM  distribution 
is  quite  irregular;  we  thus  would  not  expect  that  a  symmetric  King 
model  gives  a  good  fit. 

The  second  model  allows  for  variations  of  the  ICM  profile  with 
polar  angle.  This  model  assimilates  some  of  the  large-scale  irreg¬ 
ularities  of  the  ICM  and  allows  us  to  search  for  irregularities  on 
smaller  angular  scales.  Dividing  the  cluster  into  18  pie  sections  of 
20°  opening  angle  each,  18  King  models  are  fitted  to  the  surface 
brightness.  As  shown  in  Fig.  4  for  two  example  pie  sections,  the 
fiducial  value  of  ft  =  2/3  yields  unsatisfactory  results:  the  observed 
surface  brightness  profiles  show  a  more  pronounced  turn-down  at 
cluster  core  distances  between  5  and  10  arcmin  than  expected  for 
ft  =  2/3.  Leaving  ft  as  a  free  parameter  yields  more  satisfactory  fits. 
For  polar  angle  values  <f>  that  are  odd  integer  multiples  of  10°,  the 
surface  brightness  model  is  defined  by  the  fitted  King  profiles.  For 
intermediate  values,  we  interpolate  between  the  predicted  surface 
brightness  of  the  adjacent  two  King  models.  To  avoid  modification 
of  the  model  by  the  radio  galaxy  3C  129,  we  do  not  use  the  King 
model  of  the  corresponding  pie  section  (polar  angle  240°-260°), 
but  interpolate  between  the  King  models  of  the  two  adjacent  pie 
sections.  We  use  these  surface  brightness  models  only  to  find  ex¬ 
cess  emission;  the  statistical  significance  of  excess  features  on  small 
angular  scales  is  rigorously  evaluated  based  on  count  statistics  as 
described  above. 

The  residual  surface  brightness  distributions  smoothed  with  a 
Gaussian  of  18.8  arcsec  FWHM  are  shown  in  Fig.  3,  centre  and 
lower  panels.  The  map  derived  from  subtracting  the  symmetric  King 
model  emphasizes  the  east-west  elongation  of  the  ICM:  two  regions 


Figure  4.  Radial  0.5-5  keV  surface  brightness  distributions.  The  left  and  right  profiles  have  been  derived  from  “pie-sections”  of  the  Chandra  data  with  polar 
angles  from  80°  to  100°  and  from  240°  to  260  ",  respectively.  The  full  lines  show  fits  of  King  models  and  the  best-fitting  values  of  the  normalization  factor:  So 
(in  1 0  8  photon  cm-2  s-1  arcsec-2),  core  radius  rc  (in  arcmin),  and  /S'  are  given  in  the  upper  right  comer  of  each  panel  as  PI,  P2,  and  P3,  respectively.  In  the 
left-hand  panel,  the  dotted  line  shows  the  fitted  profile  of  the  right-hand  panel,  and  emphasizes  the  difference  of  the  two  profiles.  In  the  right-hand  panel,  the 
dotted  line  shows  a  King  model  with  ft  =  2/3. 
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with  excess  emission  can  be  recognized:  one  is  located  at  and  toward 
the  east  of  the  cluster  centre,  extending  over  at  least  7  arcmin.  The 
second  one  is  located  west/south-west  of  the  cluster  centre.  Note  that 
the  relative  strength  of  these  two  excess  regions  depends  strongly 
on  the  choice  of  the  location  at  which  we  centre  the  symmetric 
King  profile.  Moving  this  location  by  1.3  arcmin  toward  the  west 
of  the  cluster,  reduces  the  western  excess  substantially,  especially 
the  bright  region  around  the  head  of  the  radio  galaxy  3C  129.  Based 
on  this  surface  brightness  distribution,  it  is  tempting  to  speculate 
that  the  excess  region  near  the  cluster  centre  shows  remnants  of 
a  sub-cluster  passage  through  the  cluster  centre.  The  curvature  of 
the  3C  129  radio  tail  may  be  caused  by  a  sub-cluster  that  attracts 
the  radio  galaxy  and  whose  dark  matter  centre  of  gravity  even  lies 
in  the  south-west  of  the  south-western  excess  region.  Upcoming 
observations  of  3C  129  with  the  XMM-Newton  satellite  will  shed 
more  light  on  the  nature  of  these  excess  regions.  The  map  obtained 
by  subtracting  the  18-slice  pie  model  from  the  surface  brightness 
distribution  does  not  show  evidence  for  very  sharp  surface  brightness 
gradients  from  large-scale  structure  shocks  or  adjacent  ICM  phases 
of  widely  different  X-ray  emissivities. 


5  X-RAY  SURFACE  BRIGHTNESS  NEAR  THE 
TWO  RADIO  GALAXIES 

5.1  X-ray  emission  from  3C  129.1 

While  the  aimpoint  of  the  ACIS  I  observation  is  located  near  3C 
1 29. 1 ,  that  of  the  ACIS  S  observation  lies  near  the  head  of  3C  1 29  at 
about  12  arcmin  from  3C  129.1.  Thus,  although  the  latter  observa¬ 
tion  has  a  three  times  longer  integration  time,  its  angular  resolution 
is  substantially  lower:  0.84  arcsec  compared  with  9"(FWHM)  at 
the  aimpoint  of  the  observation  and  at  12  arcmin  off-axis,  respec¬ 
tively.  We  therefore  use  only  the  ACIS  I  observation  for  the  study 
of  localized  excess  emission  in  the  surrounding  of  3C  129.1. 

The  X-ray  emission  near  3C  129.1  is  shown  in  Fig.  5  together  with 
a  contour  plot  of  the  high-resolution  radio  data.  The  high-resolution 
radio  map  reveals  a  well-defined  northern  jet,  suggesting  that  the 
northern  side  approaches  the  observer.  An  excess  of  X-ray  emission 
over  that  from  the  ICM  can  be  recognized  up  to  10  arcsec  from 
the  radio  core.  The  excess  emission  corresponds  to  a  luminosity  of 
1.5  x  1041  erg  s-1.  There  is  clear  evidence  that  at  least  some  of 
this  emission  comes  from  localized  sources  rather  than  from  diffuse 
gas.  At  about  4  arcsec  to  the  north  of  the  radio  core  (source  ‘C’)  we 
detect  a  total  7  counts  (6  counts)  from  a  sky  region  of  1 . 1  arcsec  (0.4 
arcsec)  radius,  where  0.25  counts  (0.04  counts)  had  been  expected 
based  on  an  annular  comparison  region  of  1  arcsec  inner  radius  and 
20  arcsec  outer  radius.  In  Table  4  we  list  the  locations,  statistical 
chance  probabilities,  and  fluxes  for  the  extended  ‘centre  source’ 
as  well  as  for  5  potential  point  sources.  Using  the  Galactic  neutral 
hydrogen  column  density,  we  estimate  a  flux  of  6.9  x  1040  erg  s-1 
for  the  centre  source  and  fluxes  between  8.7  x  1039  and  2.1  x 
104°  erg  s-1  for  the  point  source  candidates.  The  0.5-5  keV  energy 
spectrum  of  the  combined  excess  X-ray  photons  can  be  described  by 
a  power  law  of  spectral  index  a  =  1.7+q4  (with  EdN /dE  oc  E~a). 
Using  only  the  events  from  the  point  source  candidates,  the  best¬ 
fitting  spectral  index  is  a  =  1 .3 49  4.  Comparing  the  mean  source 
spectrum  luminosity  with  that  of  resolved  sources  in  other  X-ray 
faint  galaxies  (Sarazin,  Irwin  &  Bregman  2000;  Blanton  et  al.  2001a; 
Irwin,  Sarazin  &  Bregman  2002),  the  luminosity  of  the  3C  129.1 
sources  seems  to  be  relatively  high,  while  their  energy  spectrum 
is  relatively  soft.  This  is  especially  true  for  source  C  which  has 


07.6  06.8  4:50:06.0 
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Figure  5.  The  grey-scale  plot  shows  the  0.5-5  keV  X-ray  surface  brightness 
at  the  location  of  the  radio  galaxy  3C  129.1,  smoothed  with  a  Gaussian  of 

1 . 1  arcsec  FWHM.  The  grey-scale  is  linear  and  peaks  near  the  radio  core 
of  3C  129.1  with  2.5  x  10  photon  cm-2  s-1  arcsec-2;  for  source  A  the 
surface  brightness  peaks  at  1  x  1 0  7  photon  cm-2  s-1  arcsec-2.  The  radio 
contours  have  been  derived  from  the  8-GHz  map  with  a  beamwidth  of  0. 
81  x  0.  76  arcsec2  FWHM  with  9  logarithmically  spaced  levels  from  2  to 
512  mJy  beam-1.  The  high  resolution  radio  map  reveals  a  well-defined 
northern  jet,  suggesting  that  the  northern  side  approaches  the  observer.  The 
X-ray  photon  excesses  discussed  in  the  text  have  been  marked  by  circles. 

the  highest  statistical  significance,  the  highest  flux  and  a  very  soft 
energy  spectrum  (see  Table  4). 

5.2  X-Ray  surface  brightness  modification  near  3C  129 

The  evidence  for  X-ray  emission  from  the  3C  129  radio  core  and 
jet  has  been  described  in  Paper  II.  The  excess  emission  was  well 
localized  and  came  from  a  sky  region  a  couple  of  arc-seconds  in 
diameter,  located  near  the  3C  129  radio  core.  In  Figs  3  and  6,  we 
excluded  this  core  and  jet  emission  from  the  analysis,  by  setting  the 
residual  surface  brightness  of  these  sky  regions  to  zero.  The  residual 
X-ray  map  of  Fig.  3  (lower  panel)  and  the  enlarged  view  of  the  sky 
region  near  the  radio  head  of  3C  129  in  Fig.  6  reveal  some  evidence 
for  an  excess  from  a  12  x  36  arcsec2  (6  kpc  x  17  kpc)  large  region 
at  about  30  arcsec  (14  kpc)  toward  the  east  of  the  radio  core.  Us¬ 
ing  an  annulus  of  45  arcsec  inner  radius  and  90  arcsec  outer  radius 
centred  on  the  radio  core  as  comparison  region,  we  estimate  the  sta¬ 
tistical  significance  of  the  excess  to  be  3.7 a  (ON:  131  counts,  OFF: 
92  counts). 

The  residual  surface  brightness  map  of  Fig.  3  (lower  panel)  does 
not  show  any  evidence  for  excess  or  deficit  emission  along  the  tail 
of  the  radio  galaxy  3C  129.  Indeed,  in  absence  of  strongly  emitting 
sheets  of  compressed  plasma,  we  do  not  expect  to  detect  any  excess 
or  deficit,  given  the  small  volume  of  the  tail.  Using  a  rectangular 
test  region  with  a  solid  angle  coverage  of  1.5  x  6  arcmin2  we  can 
determine  the  ICM  surface  brightness  to  a  1  a  accuracy  of  between 
4  per  cent  and  15  per  cent  near  the  radio  head  and  near  the  extreme 
end  of  the  radio  tail,  respectively.  Assuming  that  (i)  the  tail  is  fully 
depleted  of  ICM.  (ii)  it  lies  in  the  plane  of  the  sky,  (iii)  its  depth 
along  the  line  of  sight  equals  its  width  perpendicular  to  the  major 
tail  extension,  and  (iv)  that  the  ICM  can  be  described  by  a  King 
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Table  4.  X-ray  sources  of  3C  129.1. 


Name 

RA  (Eq.  2000) 
[hh:mm:ss.ss] 

Dec.  (Eq.  2000) 
[dd:mm:ss.ss] 

Non 

>]N  OFF'" 

Chance 

Probab. 

£0.5-5keV* 

[keV] 

Fluxc 

[erg  cm-2  s-1] 

Luminosity^ 
[erg  s"1] 

A 

04:50:07.  075 

45:03:02.91 

3 

0.40 

8.0  xlO-3 

2.17  ±0.29 

2.0  10-15 

8.1  x  1039 

B 

04:50:06.729 

45:03:02.92 

4 

0.40 

7.9  x  10-4 

3.58  ±  0.68 

2.7  10-15 

1.1  x  104° 

C 

04:50:06.618 

45:03:09.18 

7 

0.40 

2.3  x  10-7 

2.27  ±0.17 

5.0  10-15 

2.1  x  104° 

Centre 

04:50:06.584 

45:03:05.92 

25 

2.80 

6.7  x  10-16 

2.48  ±0.15 

1.7  10-14 

7.0  x  104° 

D 

04:50:06.149 

45:03:01.35 

5 

0.40 

6.2  x  10-5 

3.69  ±  0.58 

3.5  10-15 

1.5  x  104° 

E 

04:50:05.740 

45:03:03.61 

6 

0.40 

4.1  x  10-6 

4.16  ±0.46 

4.2  10-15 

1.8  x  104° 

a Expected  number  of  ‘background’  counts  in  the  ON  region. 
^Mean  energy  of  excess  0.5-5  keV  Chandra  photons. 
cObserved  flux. 

^ ‘Unabsorbed  luminosity’,  assumes  Nh  =  8.8  x  1021  cm-2. 
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Figure  6.  The  grey-scale  plot  shows  the  residual  0.5-5  keV  surface  bright¬ 
ness  distribution  near  the  head  of  the  radio  galaxy  3C  129  after  subtraction 
of  the  18-slice  pie  model,  smoothed  with  a  9.4  arcsec  FWHM  Gaussian.  The 
grey-scale  is  linear  and  peaks  near  the  radio  core  of  3C  129  with  a  value  of 
7.4  x  1 0  7  photon  cm-2  s-1  arcsec-2;  the  two  bins  at  the  bottom  left  side  in 
a  similar  tone  of  grey  have  a  surface  brightness  of  5.7  x  !  0  7  photon  cm-2 
s-1  arcsec-2.  The  radio  contours  have  been  derived  from  the  5-GHz  VLA 
data  with  1.8  arcsec  FWHM  beamwidth  and  the  contours  are  logarithmically 
spaced  from  1  to  32  mjy  beam-1.  An  X-ray  bright  region  southeast  of  the 
radio  core  can  be  recognized.  The  rectangle  shows  the  region  of  the  sky  that 
we  used  to  determine  the  statistical  significance  of  the  excess. 


model  with  rc  =  7.1  arcmin  and  /3  =  2/3  up  to  a  cluster  core  distance 
of  40  arcmin,  we  expect  a  reduction  of  the  ICM  surface  brightness 
of  only  between  3  and  4  per  cent  for  locations  near  the  radio  head 
and  near  the  extreme  end  of  the  radio  tail,  respectively.  Although 
rc  depends  strongly  on  the  location  chosen  as  cluster  centre,  the 
inferred  reductions  in  the  ICM  surface  brightness  do  not. 


6  PRESSURE  BALANCE  ALONG  THE  RADIO 
TAIL  OF  3C  129 


In  this  section  we  compare  the  thermal  pressure  of  the  ISM  and  ICM 
with  the  minimum  non-thermal  pressure  of  the  3C  129  radio  tail. 
The  energy  density  of  the  synchrotron  plasma  (including  magnetic 
field  energy)  is  minimized  by  the  magnetic  field  (see  Pacholczyk 
1970): 

Bmin  =  5>7t)2/7(l  +  kf'U^V-2'1  cmL s,i 


The  ^-factor  is  the  contribution  to  the  non-thermal  pressure  from  rel¬ 
ativistic  protons  over  and  above  that  from  the  relativistic  electrons, 
<f)  is  the  volume-filling  factor  for  the  radio-emitting  plasma,  and  V 
is  the  emitting  volume,  well  constrained  by  the  observations.  The 
sum  runs  over  the  considered  frequency  intervals  i  of  the  observed 
radio  spectrum,  the  value  of  Cnj  depends  weakly  on  the  synchrotron 
spectral  index  and  the  lower  and  higher  cut-offs  of  the  ith  frequency 
interval,  and  Ls>i  is  the  corresponding  synchrotron  luminosity.  We 
compute  true  minimum  energy  estimates  by  assuming  a  magnetic 
field  of  Smi„,  a  volume  filling  factor  of  (f>  =  1,  no  pressure  contri¬ 
bution  from  electrons  outside  the  energy  range  probed  by  the  radio 
data,  and  no  pressure  contribution  from  relativistic  protons  ( k  —  0). 
The  finite  angular  resolution  of  the  radio  observations  might  overes¬ 
timate  the  emission  volume  near  the  cores  of  the  two  radio  galaxies 
and  thus  further  lower  the  pressure  estimate.  The  non-thermal  par¬ 
ticle  energy  density  is  given  by 


4  B2 

T?  _  _  ^min 

p  ~  3  8tt 


(9) 


and  the  equation  for  the  minimum  isotropic  pressure  of  the  syn¬ 
chrotron  plasma  (including  magnetic  field  pressure)  reads 


1  9 

p  .  —  -p  j _ i m 

*  min  —  0  i  n  o 

3  3  87t 


(10) 


In  the  latter  equation  we  conservatively  assumed  a  tangled  magnetic 
field. 

Based  on  these  standard  assumptions,  Pmin  can  be  computed  from 
the  radio  data  alone.  As  a  starting  point  we  use  either  the  0.33-,  5- 
or  8-GHz  surface  brightness  radio  map.  For  each  point  along  the 
radio  tail  we  compute  using  thin  cylindrical  segments,  in  other 
words,  discs,  seen  edge  on.  The  disc  axes  are  aligned  with  the  radio 
tail.  The  width  (and  depth)  of  a  disc  is  given  by  the  two  points  at 
which  the  radio  intensity  drops  to  10  per  cent  of  its  peak  value. 
For  each  disc  we  compute  the  mean  radio  flux  So  at  the  frequency 
Vo  of  the  considered  map,  by  averaging  the  surface  brightness  over 
the  cross-section  area  of  the  disc.  Given  So,  we  construct  the  0.33- 
8  GHz  radio  energy  spectrum,  assuming  the  measured  0.33-1 .4  GHz 
and  5-8  GHz  spectral  indices  over  the  frequency  ranges  from  0.33- 
2.65  GHz  and  2.65-8  GHz,  respectively.  We  use  the  ICM  pressure 
front  the  spectroscopic  ICM  deprojection  analysis  of  the  Chandra 
data  (Paper  I)  and  neglect  the  pressure  contribution  of  the  ICM 
magnetic  field.  Taylor  et  al.  (2001)  estimate  that  the  ICM  magnetic 
field  might  be  as  strong  as  10  pG.  If  this  is  true,  the  thermal  and 
magnetic  ICM  pressures  would  be  comparable  at  and  near  the  head 
of  the  radio  galaxy  3C  129. 
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Figure  7.  The  thermal  ICM  and  non-thermal  radio  plasma  pressures  for 
the  radio  galaxy  3C  129  as  function  of  distance  to  the  cluster  centre  at 
a?  =  04h50m02l9  and  S  =  45°02'26'.'0  (J2000).  The  thermal  ICM  pressure 
from  the  spectroscopic  deprojection  analysis  of  the  X-ray  data  (Paper  I)  is 
shown  by  the  solid  line  and  the  data  points.  The  dashed  and  dotted  lines  show 
the  minimum  pressure  of  the  radio  plasma  for  the  northern  and  southern  radio 
tail,  respectively,  computed  for  a  minimum  radio  frequency  ui  =  330  MHz. 
An  estimate  of  the  pressure  contribution  from  very  low  energy  electrons 
is  shown  by  the  dash-dotted  lines,  which  have  been  computed  under  the 
assumption  that  the  radio  spectrum  extends  down  to  Vi  =  1  MHz  with  the 
same  spectral  index  as  measured  for  330  MHz-1.4  GHz. 

The  pressure  profiles  are  shown  in  Fig.  7.  Over  the  first  13  arc- 
sec  (6.1  kpc)  from  the  nucleus  of  3C  129  the  minimum  pressure 
of  the  radio  plasma  exceeds  the  pressure  of  the  ICM  by  a  factor 
of  up  to  ~5,  suggesting  either  self-confinement  or  confinement  by 
the  cold  disc-like  remnant  of  the  ISM  that  has  been  stripped  off  the 
galaxy  by  the  relative  motion  of  galaxy  and  cluster  gas  (Toniazzo 
&  Schindler  2001).  At  2  arcmin  (57  kpc)  from  the  radio  core,  the 
minimum  pressure  fell  by  more  than  one  order  of  magnitude  and 
shows  strong  wiggles,  suggesting  ‘stormy  local  ICM  weather’.  Un¬ 
fortunately,  low  X-ray  photon  statistics  do  not  allow  us  to  determine 
the  ICM  properties  with  sufficient  statistical  accuracy  to  assess  the 
imprint  of  this  stormy  weather  in  the  X-ray  surface  brightness  or 
X-ray  energy  spectra.  Further  away  than  2  arcmin  front  the  radio 
core,  the  minimum  pressure  levels  off  at  roughly  l/40th  of  the  ICM 
pressure. 

The  equipartition  magnetic  fields  lie  between  16  and  1  pG  near 
the  beginning  and  the  end  of  the  tail,  respectively.  The  lower  fre¬ 
quency  cut-off  of  330  MHz  corresponds  to  Lorentz  factors  of  y\— 
4000(6/5pG)_I/2.  We  estimate  the  possible  pressure  contribution 
of  electrons  with  very  low  Lorentz  factors  by  extrapolating  the  radio 
energy  spectrum  down  to  iq  =  1  MHz,  corresponding  to  Lorentz 
factors  of  y  i  ~  220(6/5  p.GR1/2.  Using  the  lower  frequency  cut¬ 
off,  the  inferred  equipartition  magnetic  field  strength  increases  to 
between  22  and  3  pG  near  the  beginning  and  the  end  of  the  tail, 
respectively.  The  corresponding  pressure  estimates  are  shown  by 
the  dashed  dotted  lines  in  Fig.  7. 

For  the  first  10  arcmin  of  the  tail,  the  radio  spectrum  is  suffi¬ 
ciently  flat  that  the  lower  cut-off  frequency  does  not  qualitatively 
alter  the  result:  the  minimum  non-thermal  pressure  lies  substantially 
below  the  ICM  pressure.  We  will  discuss  the  implications  in  the  next 
section. 


7  SUMMARY  AND  DISCUSSION 

The  Chandra  observations  show  that  the  radio  galaxy  3C  129.1  is 
associated  with  an  extended  X-ray  excess.  The  excess  can  be  di¬ 
vided  into  5  point  sources  and  a  central  component.  On  the  small 
angular  scale  of  the  central  component  (4  arcsec  diameter)  the  ISM 
is  not  expected  to  contribute  to  the  excess.  We  therefore  consider  it 
most  probable  that  the  central  excess  emission  is  produced  by  galac¬ 
tic  sources  and  possibly  also  by  the  central  active  galactic  nucleus 
(AGN). 

For  the  radio  galaxy  3C  129  we  find,  in  addition  to  the  core  and 
jet  emission  described  in  Paper  II,  some  evidence  for  extended  X- 
ray  emission  from  a  12  x  36  arcsec2  (6  kpc  x  17  kpc)  large  region 
to  the  east  of  the  radio  core.  If  confirmed  with  higher  statistical 
significance,  a  natural  explanation  of  the  excess  would  be  thermal 
emission  from  plasma  compressed  and  heated  by  the  bow  shock  that 
precedes  the  galaxy. 

Comparing  the  minimum  pressure  of  the  3C  129  radio  tail  with 
the  thermal  ICM  pressure,  the  radio  plasma  appears  to  be  under¬ 
pressured  by  a  factor  of  between  10  and  40.  Over  the  first  10  arcmin 
of  the  radio  tail,  the  radio  spectrum  is  sufficiently  flat  so  that  the 
results  depend  only  weakly  on  the  assumed  minimum  Lorentz  factor 
of  the  non-thermal  electron  distribution. 

Remarkably,  a  very  similar  pressure  mismatch  has  been  found  for 
the  kpc  radio  features  of  a  number  of  FR I  and  FR  II  radio  galaxies 
(see  Hardcastle  &  Worral  2000,  and  references  therein). 

For  most  radio  galaxies  the  magnetic  field  inside  the  radio  plasma 
is  poorly  constrained  toward  higher  values  and  might  dominate  the 
plasma  pressure.  In  the  case  of  3C  1 29,  we  have  an  additional  handle 
on  the  magnetic  field  strength  from  the  constraints  on  the  radiative 
age  of  the  radio  electrons:  If  the  missing  pressure  of  the  radio  plasma 
is  provided  by  the  magnetic  field  alone,  the  implied  field  strengths 
of  between  15  pG  and  30  pG  translate  into  very  high  lower  limits 
on  the  galaxy  velocity  of  between  1 1  000  km  s-1  and  31  000  km  s-1 
(see  equation  6).  Buoyant  motion  of  the  radio  plasma  would  proceed 
with  a  velocity  lower  or  equal  to  the  ICM  sound  speed,  and  would 
thus  not  qualitatively  alter  our  conclusions.  The  apparent  pressure 
mismatch  thus  suggests  one  of  the  following  possibilities:  (i)  ther¬ 
mal  or  relativistic  protons  and/or  thermal  or  relativistic  low-energy 
electrons  carry  the  major  fraction  of  the  radio  plasma's  pressure; 
(ii)  the  volume  filling  factor  of  the  radio  plasma  cj>  is  substantially 
smaller  than  unity:  $  =  0.0016-0.018;  (iii)  the  6  field  is  much  larger 
than  the  equipartition  magnetic  field  and  electrons  are  re-accelerated 
inside  the  radio  tail;  (iv)  the  radio  galaxy  lies  behind  or  in  front  of 
the  plane  of  the  cluster  so  that  the  true  ambient  ICM  pressure  is 
lower  than  assumed  in  our  analysis.  The  detection  of  an  ‘ICM  cav¬ 
ity"  associated  with  the  radio  tail  would  make  it  possible  to  estimate 
the  volume  filling  factor.  We  hope  to  detect  such  a  cavity  with  an 
upcoming  XMM-Newton  observation  of  3C  129. 
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